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Abstract 

We generalise the equations governing relativistic fluid dynamics given by Ehlers and Ellis for 
general relativity, and by Maartens and Taylor for quadratic theories, to generalised f(R) theories of 
gravity. In view of the usefulness of this alternative framework to general relativity, its generalisation 
can be of potential importance for deriving analogous results to those obtained in general relativity. 
We generalise, as an example, the results of Maartens and Taylor to show that within the framework 
of general f(R) theories, a perfect fluid spacetime with vanishing vorticity, shear and acceleration 
is Friedmann-Lemaltre-Robertson-Walker only if the fluid has in addition a barotropic equation of 
state. It then follows that the Ehlers-Geren-Sachs theorem and its "almost" extension also hold for 
f(R) theories of gravity. 
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1 Introduction 

Recently, Maartens and Taylor examined the kinematical and dynamical properties of 4-D fluid 
spacetimes in quadratic gravity — with the gravitational Lagrangian density of the type 

£ G = ^;V=g[R + aR 2 + pR^R^] , (1) 

where a and (3 are positive coupling constants, R is the Ricci curvature scalar and i? M „ is the Ricci 
curvature tensor. They derived the generalisations of the equations governing relativistic fluid dynamics 
given by Ehlers || and Ellis [^0) for general relativity (GR) to the quadratic case. The kinematic and 
dynamic formulation of fluid spacetimes has been extremely useful in studying a variety of problems 
within the context of GR and it is similarly expected to be of value in the case of higher-order theories 
of gravity, as it gives an alternative framework to study problems analogous to those considered in GR. 
An interesting example of this was given by Maartens and Taylor, when they showed that the usual 
kinematic characterisation of the Fricdmann -Lcmaitrc-Robertson- Walker (FLRW) geometry no longer 
holds in the quadratic theory, but that it can be restored, if in addition the fluid is assumed to be 
barotropic. A consequence of this is that the Ehlers-Geren-Sachs theorem (EGS) || also holds in the 
setting of quadratic theories of gravity. The EGS theorem states that if a family of freely falling observers 
measure the self-gravitating background radiation to be everywhere exactly isotropic, then the Universe 
is exactly homogeneous. The importance of this theorem lies in the fact that the observable isotropy of 
the cosmic microwave background radiation (CMBR) may be employed^ to deduce that the spacetime is 
also homogeneous and therefore of the FLRW type. 

A shortcoming of the above work is that, as is customary, it confines itself to the quadratic theories, 
thereby implicitly assuming the quadratic theories to be the generic prototypes of the higher-order theories 
in general. Higher-order corrections to the action arise in attempts at generalising GR as for example 
in (i) the low energy limit of superstrings |l , (ii) in perturbation expansions of GR, where the existence 
of such higher-order terms can eliminate the divergences that would otherwise arise ^| and (iii) in 
the renormalisation of higher-loop contributions which gives rise to terms in the action that are higher 
than quadratic. Such terms are thought to have been dominant at early epochs in the evolution of the 
Universe, when the curvature of spacetime was high. 

This indicates that the genericity of the quadratic theories cannot be assumed a priori and higher 
than quadratic terms in the action are also likely to be involved in generalisations of GR. Further, there 
is no a priori reason why results obtained within the context of quadratic theories remain stable (in 
the sense of 0) with respect to cubic (and higher-order) perturbations to the action. In fact there are 
examples where the inclusion of cubic terms to the action can produce qualitatively important changes 
to the corresponding results obtained within the framework of the quadratic theories (see jll| for an 
example). To have a feel for this, recall that general higher-order theories with gravitational Lagrangian 
density of the form 

^G = ^V=gf(R), (2) 

where f(R) is a differentiable function of R and df(R)/dR ^ 0, can be expressed as GR plus a scalar field 
<fi, with a corresponding scalar potential V(<p) || [uj. Now it turns out that the asymptotic behaviour 
of this potential depends crucially on the form of the action. For example, in 4-D spacetimes, the scalar 
potential V(<p) corresponding to quadratic theories tends asymptotically towards a finite constant limiting 
value, whereas for cubic and higher-order actions, it tends to zero Higher than quadratic terms in 
the action could therefore produce important consequences. 

As a result, it is of importance to generalise the equations governing relativistic fluid dynamics given 
by Maartens and Taylor to the case of general higher-order gravity with gravitational Lagrangian den- 
sity given by (||). This would be of value, as it allows an alternative framework within which general 
higher-order theories can be studied. Here as an example, we study the conditions for the kinematic char- 
acterisation of the FLRW models (and the EGS theorem) , to see whether they are stable with respect to 
this action. We then use our results to test the stability of the almost-EGS theorem of Jl5|, |l7j| . 

The plan of the paper is as follows. In section 2 we generalise the results of [jl3|, by deriving the 
Ehlers-Ellis equations of relativistic fluid dynamics for a general action with a gravitational Lagrangian 

1 Assuming the isotropy to hold everywhere. 
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density of the form (|2|) and a general energy- momentum tensor T pv . In section 3 we consider the FLRW 
models, the EGS and almost-EGS theorems within the context of these theories, and finally section 4 
contains the conclusion. 



2 Fluid Kinematics and Dynamics of f(R) Theories 

In this section we generalise the equations of relativistic fluid dynamics given in ||, to the case of 
generalised theories of gravity given by the action 



S = J- [ J~g f(R) d 4 x + [ C M d 4 x , (3) 



2k q 



where f(R) is some arbitrary differentiable function of the Ricci curvature scalar, Cm is the Lagrangian 
density for the matter sources and we have chosen units such that c = 1 and kq = 8nG. Throughout, 
Greek indices run from /i, v = . . . 3 and in a coordinate system {a; , x a } (a = 1, 2, 3), the spatial coor- 
dinates are denoted by Latin indices a, 6, . . .. 

The variation of this action with respect to the metric tensor of spacetime, g pv , results in the gener- 
alised field equations, given by || 

Rfj,u = f'- 1 [ \fg»u + I" (V M V„i? - g pv OR) 

+ f" (V P R V U R - V P R V P R) + k T piJ } , (4) 

where /' := df(R)/dR, etc, OR := g^V^V^R, and the energy-momentum tensor of the matter sources 
is defined by the variational derivative T» v := {2/*/=g) 5C M /5g pi/ . (Note that V p T A " / = still holds 

!■) 

The contraction of this equation gives the Ricci curvature scalar in the form 

R = f'- 1 [2f-3f"nR-3f'"V p R\7 p R + K T] , (5) 
where T denotes the trace of T pu . Taking the covariant derivative of the Ricci scalar (|^) one obtains 



M f - f" R-3f"'aR-3f""\7 p R\7PR 



Now following Ehlers [|8| and Ellis |l0[, we define a projection tensor, which allows geometrical quan- 
tities to be projected into the instantaneous rest spaces of observers co-moving with the timelike matter 
flow with (normalised) 4-velocity u p , given by h p v := 5 P V + u p u v . Projecting parallel and orthogonal to 
the matter flow lines one obtains a 1 + 3 decomposition of the dynamic equations ("threading picture"). 
In particular this yields a decomposition of the covariant derivative of the matter 4-velocity u p , which is 
used to define kinematical quantities characteristic of the matter flow, in the form 

VvUp := -a p u v + <7( MI /j + \<dh l _ ll/ + w^j , (7) 

where a p denotes the acceleration vector, (J( pv ) the rate of shear tensor, Q the rate of expansion scalar, 
and uj^] the rate of vorticity tensor. Also in line with J8|, [T^] , one can define the usual vorticity vector 
in the form (e i23 = ^ ■= \ <^ vpa u vp u a io pv = t pvpa uo p u° ' . 

The "electric" and "magnetic parts" of the trace-free Weyl conformal curvature tensor can be defined 

by 

E( plJ ){u) := C ppvc ,u p u a , H^(u) := \ e ppa p C a& v< , u p u a . (8) 



2 FLUID KINEMATICS AND DYNAMICS OF F(R) THEORIES 



Splitting the field equations as in the case of GR (and the quadratic case [|l3|), by contracting the 
Ricci tensor (0) with the 4- velocity u p and the projection tensor h%, yields 



f _1 [ -| / + /" W v V M V V R + f" h>» V P R V V R + k T pv u» u v ] 



Rv P h\u p = f'- 1 yf"h\{V v V p R)u p + f"'Rh\V v R + K T vp h\u p 
RpahWv = f- 1 [\fh,„, + f"h^h a u(y p V a R-h pa uR) 

+ f" h% h% (V P R V a R - V V*i? V A i?) + k T pa h p p h% ] . 



(9) 
(10) 

(11) 



As in the case of GR || ^0| , starting from the Ricci identities for the vector field u p , multiplying by u v 
and projecting into the local 3-space, gives the propagation equation for v pu = h a p hr u V Q it/3, in the form 



h a p h,P v v a p - a p a v - h a p h^ v Va^ + v p p v vp + R ai ps h a p u 1 h l} u u s = . 
The trace of this equation yields the generalised Raychaudhuri equation for the action ([|) 



e 
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9^ - V p a p + 2 (o- 2 - lo 2 



+ [ -y + f"h^V p V„R + f'"h^V p RV„R + K oT^u p u» ] . 
The skew part of (jlj) is the generalised vorticity propagation equation for the f(R) theories, 

= h\ h 13 ^ (uj a/3 + V a ap) - 2 a p [p u v ] p + § 9 uj p „ , 



(12) 



(13) 



(14) 



which is identical to the Einstein case. On the other hand the symmetric, trace-free part of ([12]) is the 
shear propagation equation for the action (|^) given by 



= h a ( p h? v) (b a p, -V a ap) - a p a u + lu^uj^ + a pp a p v + §9 a pv 

-If'- 1 (h% - \ V h at3 ) [ /" W a W p R + f" V a RV p R + n T aP } 



(15) 



One can also obtain the three constraint equations from the Ricci identities for u p in the form 
= ^(fVe-^V/') + e pvpa u a (V l/ u Jp -2u Jv a p ) 



which are the field equations fllO|), together with 



/" hP v (V U V P R) u p + f" R h p v V U R + K T% h% u p 



h^v V p uj L/ = a p u p 
H pv = 2 a {p uj v) + h% h p v) (VW + V 7 <r Q 4 ) e^ S e 



(16) 

(17) 
(18) 



Using the decomposition of the Riemann curvature tensor into the trace- free Weyl conformal curvature 
tensor, the Ricci curvature tensor and the Ricci curvature scalar, 

R P v P a • — C pvpa -\- 9 P \ P R(t\v ~t~ 9v[<7 3 ^ 9jj\p9(r\v : (19) 

the Bianchi identities 

VfrRvpiar = (20) 

can be cast into the form |12] ] 

VaC ^P<? _ _ v ln R v]p _ i gPlnyv] R . (21) 
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Now using equations (|6|)-(|ll|), the constraint and evolution equations for Ep V and Hp V can be derived 
for the case of f(R) theories, in the form 



^ a01 o a s H 6 j3 u 7 - 3H" a u a 



R R a p u p - (R Pl u + \R)V a R 
V Q (DR) - (VVpRYuP + hl\ V" (VpV^R) ] , 



(22) 



h" a Wp E a0 



3EK a v)a + EK u")" -QEP v - U LV E a p a a/3 
- 2HK e v)a ^ ap u 7 + hK e v ^ s V pH^ u s 



K fl—1 



+ if/ 



hK h% v°r" 7 ui - hK h% (T ap ) ■+\fh tiv 



h" a h~<p V 7 ff a/3 



+ \f'- 1 f"h^ a h% [r{R^-\R 9 
+ V" (V 7 V /3 i?) u 1 - (V Q V P R) ' ] 
- e^ 7 a a5 E 5 p u 7 + 3E" a uj a 



ap\ 



«0 jl -1 e ^a(3~{ 



V a T p u 7 US 



h» a h"pH 



ap 



+ i /' -i f ^fh [ Vq V^V 5 i? - V Q i? ^ }u 7 u s , 
3HK cr v)a + HK v v)a - 6 H» v - H aP a a(3 
+ 2EK f? )afh a p u 7 - hK t v)M8 VpE a 7 u s 
+ ff'^hK^^VpT^us 



(23) 



(24) 



+ U'- 1 f"hK 



M)p-f& 



[ VpRR% - V /3 V 7 V a i? ] u 5 . 



(25) 



These equations generalise the analogous equations given for the quadratic Lagrangian theories of gravity 
fl3f to the f(R) case with a general energy- momentum tensor T^„. 

And finally we derive, by use of the Gaufi equation, which holds for spacetime configurations with 
a hypersurface-orthogonal matter fluid flow (u^„ = 0), the 3-Ricci curvature tensor of the spacelike 
3-surfaces orthogonal to w M . The 3-Riemann curvature tensor 3 Rp Vpa is given by 



Rnvpa -RpP R/jrr h h 



Rftr hfJ.p Rvp h>[M(T 2 ^ (^PP huo 



hyp} 



Using the field equations (fll]), we can split the 3-Ricci tensor 3 Rp V into its trace 



''R 



3 w 



/ + /" (2 h^VfrfuR - 3 OR) + /"' (2 R 2 - V^V-R) + n a (2 T^vPu" + T) 



2 a 2 

and its trace-free part 

3 Rp V - \ 3 Rhp V = h\h v )(V ' a ap - & a p) - Q<7p V + a M a„ - \hp V XI p a p 



(26) 
(27) 



(28) 



+ /' - 1 (h% hK - § V h aP ) [ f" V a VpR + f" V a R VpR + k T a p ] 



These relations then complete the generalisations of the kinematical and the dynamical equations given 
by Ehlers and Ellis g for GR and those given by |D| for the quadratic gravity to the case of f(R) 
theories with a general energy- momentum tensor Tp V . 



3 FLRW Cosmological Models in f(R) Theories 

In this section we look for conditions for the kinematic characterisation of FLRW geometry in the context 
of f{R) theories. To begin with let us recall that the kinematic conditions, in the context of GR and for 
perfect fluid spacetimes, are given by 



fl p — W(U/ — CT py — , 



(29) 
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which with coordinates x a — (t,x a ) allow the metric to be expressible in the following form 

ds 2 = -dt 2 + A(t) 2 A ab (x c ) dx a dx b , (30) 

such that the spacelike 3-surfaces defined by the metric A ab (x c ) have constant curvature. In usual 
Einstein's general relativity, the conditions (£9j) along with the constraint equations (|l6|) - ( |l8| ) imply 
the vanishing of the Weyl tensor. However, for the f(R) generalised equations with perfect fluid matter 
source, which we shall assume here, only is immediately zero, and therefore to show that E^ = we 
need to show that the Ricci scalar is spatially homogeneous. To do this, we recall Jl]] that the kinematic 
conditions ( 29 ) imply the existence of the metric 

ds 2 = -dt 2 + S(t,x c ) 2 \ ab (x c )dx a dx b it" =£"(), (31) 

where 5 © = VoS/S 1 . The (0a) field equations (|ic| ) together with the constraint equations (16) imply 
that 

d a do (In S) = § V a 6 = - V a V (/') . (32) 

Now, if we assume that our matter fluid source obeys a barotropic equation of state of the form p — 
p(p), dp(p)/dp 7^ 0, where p and p denote the isotropic pressure and the total energy density of the 
matter fluid respectively, then the conservation equations V ' U T^ U — with the kinematic conditions ( p9| ) 
imply that 

h\ V uP = h v ^ vP = h v ^ v Q = . (33) 

In the coordinates of (|3l|), ( |33|) gives = 0(t), therefore, ( |32"|) implies S = A(t) B(x a ), which reduces 
the metric @ to the form ©. 

The next step is to show that A ab (x c ) has constant curvature. Again, from ( |32j ) and (^3|) we have 
V a Vo (/') = which, using the connection derived from (|30|), integrates to 

f = A(t)g(x a )+B(t) , (34) 

where g and B are arbitrary functions of x a and t respectively. To proceed from here is not as straightfor- 
ward as the proof in We use equations (5) and (9) to eliminate the term involving /. The resulting 
equation, using (30), (31) and the perfect fluid form of the energy-momentum tensor is 



6A 2 6A ''/,'• 1 . ..,„_ _ .., _ .., 

(35) 



A 2 A A 2 f 



-1 A ab V a V b (/') + 3 V Vo (/') + 3k (p + p) 



where 3 R*(x c ) is the 3-Ricci scalar of the spacelike 3-surfaces defined by the metric A a b(x c ). Using 
equation (34) and the form of the metric (30), equation (35) can be reduced to (provided A ^ 0) 

h(x a ) + B(t) 3 R*{x a ) + 9A(t) 2 g(x a ) + C(t) = (36) 

where h(x a ) and C(t) are given functions and A, B and g are defined by (30) and (34). It is clear from 
equation (36) that V a g = and consequently that V a (/') = 0, which in turn implies 

h%V v R = . (37) 

Now to show that A a b(x c ) has constant curvature we proceed as follows. Equation (|37]) and the (GR) 
FLRW kinematic conditions (^9|) imply [|l] 

V^V^i? = Ru^ u u — (hA/Ajh^, (38) 

which by (15) and (28) implies 

Epv = = 3 R l _ tu — 3 3 Rhf JiV , (39) 

where 3 R is the 3-Ricci tensor formed from h ao (t, x c ), proving that the spacelike 3-surfaces t = constant 
are intrinsically isotropic. And finally we need to show that such 3-surfaces are isotropically embedded. 
The reasoning for this remains unchanged to the one given in [|l3| for the quadratic theories. First of all 
recall that Q 

R ab = 3 R\ b {x c ) + (AA + 2A 2 ) A ab , (40) 
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where 3 R* a b(x c ) is the 3-Ricci tensor formed from A t,(x c ). It can also be shown, using the field equation 
(11) together with (9) and (5), that 

R ab = (\A 2 R~AA) K ab . (41) 

Now since A and R are functions of t only, equations ( [Io| ) and (fill ) imply that 3 R* a b(x c ) — 2K A a b(x c ), 
where 

2K=\A 2 R-2AA~2A 2 = constant , (42) 

and K is the constant curvature of A a b(x c ). This then shows that the spacelike 3-surfaces t — constant 
are isotropic and isotropically embedded and therefore the metric is FLRW. The corresponding Bianchi 
identities ( p2[ ) - ( p5| ) are all identically satisfied, except for the £^„ propagation equation ((2^), which 
reduces to the field equation (11). 

In this way we have a generalisation of the result given in p3[ , thus: 

In the general f(R) theories of gravity, a perfect fluid spacetime with vanishing vorticity, shear and ac- 
celeration is FLRW only if the fluid has in addition a barotropic equation of state. 

An important consequence of this result is that the Ehlers-Geren-Sachs theorem also holds for the 
f{R) theories of gravity. 

In this sense then, the main conclusions of Maartens and Taylor remain stable to any perturbations 
of the action. This represents a significant factor in favour of the higher-order theories given that the 
EGS theorem is one of the key results motivating the standard big-bang model of the Universe. 

Finally, it has recently been shown [fl5| that the observed almost-isotropy of the cosmic microwave 
background radiation (CMBR) implies that the Universe has been almost spatially homogeneous and 
isotropic since the photon decoupling In terms of covariant perturbation theory fl5fl , this means that 
deviations of the spacetime kinematics and dynamics from FLRW are at most of order one with respect 
to some smallness parameter e. In other words, if a quantity disappears in a FLRW spacetime, it is of 
order one (or higher) in an almost-FLRW spacetime, which allows the metric to be written in a perturbed 
FLRW form Subsequently, Taylor and Maartens JlTj have shown that this theorem also holds in 

quadratic theories of gravity and may have applications to the very early Universe. 

Now the results of the section 2, together with the above result that EGS theorem also holds for the 
general f(R) theories of gravity, imply that the almost-EGS theorem shown to hold for quadratic gravity 
of |lj is further generalisable to the case of f(R) theories. We see this as follows. 

We use the formalism of |ll], [TtJ and elaborate only on the deviations from the proof of jl7) where 
necessary. We consider an acceleration- free (a M = 0) expanding congruence (8 > 0), where the kinematics 
is determined by the shear, vorticity and expansion. It can be shown in a manner analogous to that in 
flEf , that the components of the energy-momentum tensor, which deviate from perfect fluid form, are not 
dissipative quantities. Then the energy flux and the anisotropic stress, and ir^, measure the deviation 
of the distribution function from isotropy and, since they vanish in FLRW spacetimes, are first order 
quantities, i.e. 

<Z M , *m» = 0[1] ■ ( 43 ) 
From the conservation equations it is possible to derive the kinematic conditions ( p9| ) to first order, 
giving 

= 0[1] . (44) 



Furthermore, it is also possible to establish the conditions ( |33| ) to first order: 

h\ V„ P , h^Vup , h^VuQ = 0\l] . (45) 



Equations (|43j), fl44| ) and ([15| ) arc crucial to proving almost-homogeneity of the scalar curvature. From 
the constraint equations ( |16[ ) we have that 

(VV) u " = 0[l] ■ (46) 



"This amounts to an almost EGS theorem. 
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By taking the spatial derivative of the generalised Raychaudhuri equation ( |l3| ) and using (|46| ) we obtain 

/"(8 9 + 2e 2 -2i?) -/' h" M V„R = 0[1] , (47) 

which in view of the fact that 7?, Q and are non-zero in the exactly isotropic (FLRW) case, they would 
be O[0] and we can deduce 

h^VvR = 0[1] . (48) 

Then from |[|) and @ we obtain 

, = 0[1], (49) 

which establishes the almost-FLRW kinematics to 0[1]. Equation (gq) implies also that the 3-surfaces 
of constant energy density are almost-isotropic and almost-isotropically embedded. We see this from 
equations (28) and (42), where it follows from ( Eg ) that 

3 R pv - | 3 Rh^ = 0[l] (50) 

and 

2K= \A 2 R-2AA-2A 2 , K = 0[1\. (51) 

Thus, it is possible to recover all the standard relations governing an FLRW model of the Universe to 
O[0], and in particular, to show that there exists an almost-FLRW metric in general f(R) gravity. This 
then shows that the "almost" EGS generalises to the case of f(R) theories. 



4 Conclusion 

We have derived the generalised kinematic and dynamic equations governing relativistic fluid dynamics 
(given by Ehlers and Ellis for GR and for the quadratic theories by Maartens and Taylor) in the case of 
generalised f(R) theories of gravity and with a general T^ v . Such generalisation is of potential importance 
as it gives an alternative framework in order to derive results analogous to those obtained in GR for the 
generalised theories of gravity. 

As an example of application of these results, we have generalised the results of Maartens and Taylor 
to show that within the framework of f(R) theories a perfect fluid spacetime with vanishing vorticity, 
shear and acceleration is FLRW only if the fluid has in addition a barotropic equation of state. It then 
follows that the Ehlers-Geren-Sachs theorem holds for general f(R) theories as well. In this sense, then, 
the results of Maartens and Taylor are stable to perturbations of any order in the action. We also show 
that the "almost" EGS theorem of Stoeger, Maartens and Ellis and its generalisation derived by Taylor 
and Maartens in the case of quadratic gravity also remains stable to perturbations to the Lagrangian of 
the form given in equation (3). 
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